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IMPROVEMENTS TO THE SADSAM COMPUTER PROGRAM FOR AEROELASTICITY ANALYSIS 

INTRODUCTION 

This report contains a summary description of the work accomplished under 
NASA contract NAS1-12017. In general the work includes enhancements in the 
capability and efficiency of the SADSAM computer program at LRC and illustrative 
analyses demonstrating gust analysis capabilities of SADSAM. 

The modifications to the SADSAM computer program included: 

1. Incorporation of 3 optional methods of normalization for normal mode 
shapes and complex eigenvectors. 

2. The addition of the capability to determine the element and aerodynamic 
forces from complex eigenvectors. 

3. The addition of the capability to determine the element and aerodynamic 
energy distribution from complex eigenvectors. 

4. The addition of the capability to determine the modal composition of 
of complex eigenvectors. 

5. The addition of the capability to determine zero frequency (rigid body) 
normal mode shapes. 

6. Enhancement of the aerodynamic formulation by the addition of provisions 
for control surfaces, wing downwash effects on tail surfaces, chord- 
wise displacement and chordwise forces. 

7. The improvement of the efficiency of transient analysis. 

8. Miscellaneous corrections and convenience improvements. 

Appendices attached to this report were prepared for the "Theoretical 
Basis of the SADSAM Computer Program." The numbering of these appendices is 
consistent with that document and therefore results in numbering gaps in this 
report. 

*The contract research effort which has lead to the results in this report 
was financially supported by USAAMRDL (Langley Directorate) . 



EIGENVECTOR NORMALIZATION 


Prior to the SADSAM program modification under this contract, normal mode 
shapes were normalized to a unit generalized mass and complex eigenvectors 
were normalized to unity at the maximum displacement magnitude. 

At the user's option complex eigenvectors also can now be normalized to 
an effective unit generalized mass as defined by equation D-2 of appendix D. 
Provision has been made, through the use of the 0MITGENM data card, to omit 
selected degrees of freedom from the calculation of the "generalized mass" 
of complex eigenvectors. This provision is made so that non-structural masses, 
such as those associated with lag functions, may be excluded from the calcula- 
tion of "generalized mass." 

Also at the user's option either normal mode shapes or complex eigenvec- 
tors can be normalized to a unit value at the degree of freedom having the max- 
imum displacement magnitude or a user selected degree of freedom. 

COMPLEX EIGENVECTOR FORCES § ENERGY DISTRIBUTION 

The SADSAM program has had the capability of calculating the elastic 
element forces and strain energy distribution for normal mode shapes. The pro- 
gram modifications made under this effort extend this capability to complex 
eigenvectors. 

The definition of force and energy distribution of a complex eigenvector 
is more comprehensive than in the real case since it includes the consider- 
ation of mass elements, damper elements, and aerodynamic forces as well as the 
elastic elements included for the normal mode shapes. Appendix D contains a deri- 
vation and discussion of the complex force and "complex energy" calculated for 
complex eigenvectors. 
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MODAL COMPOSITION OF EIGENVECTORS 


Flutter analyses are made by SADSAM in terms of the physical coordinates 
of the system. In many cases insight into a flutter phenomenon can be gained 
through the examination of the normal mode components of the complex eigenvec- 
tor associated with the flutter root. For these purposes the capability of 
calculating modal components of complex eigenvectors was added to the SADSAM 
program. 

A derivation and discussion of the method used is contained in appendix 

D. 

DETERMINATION OF RIGID BODY MODES 

Analyses of free systems which included rigid body modes has always 
been within the capability of SADSAM. In the determination of normal modes of 
free systems, only the finite frequency modes were calculated. In order to 
fully exploit the new capability to determine modal composition, it was 
necessary to add the ability to calculate rigid body mode shapes. A derivation 
and discussion of the method used is contained in appendix F . 

' ‘ ‘ REVISED AERODYNAMIC FORCE REPRESENTATION 

One of the development tasks has been to extend SADSAM’ s existing 
capability for representing aerodynamic forces on a swept -back lifting 
surface. The existing capability is described in Appendix C of the 
SADSAM Theoretical Manual. In the existing capability, the lifting 
surface is treated as a structural beam with out-of -plane bending and 
torsional degrees of freedom. The aerodynamic forces include lift and 
moment about a structural reference axis. They are computed according 
to a modified strip theory which permits the user to specify the center 
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of pressure, the lift curve shape, and the lag function parameters for 
each strip. He can also specify a gust function. The calculations are 
performed in either the time domain or in the frequency domain at the 
user’s option. Flutter analysis is performed by the p-method, using 
SAD SAM’s standard complex eigenvalue extraction technique. This method 
has the advantages that the computed decay rates are realistic, and that 
results are obtained for each specific speed without recourse to inter- 
polation, as would be required by the more commonly used reduced frequency 
method. 

The extensions which have been implemented include: 

a. Addition of degrees of freedom to represent, at each strip, a 
sealed flap without aerodynamic overhang. 

b. Inclusion of chordwise aerodynamic forces and a chordwise 
degree of freedom for each strip. The intended purpose of this 
extension is to permit the analysis of low frequency rigid body 
modes (e.g., phugoid motion) as well as the analysis of lifting 
surfaces with low chordwise rigidity. 

c. Representation of wing-induced downwash on tail surfaces (or 
canard- induced downwash on wing and tail surfaces), including 

time delay due to convection. 

The theory for each of these effects will be described separately. 

Flap Aerodynamics 

Unsteady aerodynamic forces for an airfoil with flap in an incom- 
pressible, two-dimensional flow are derived in Theodorsen’s classical 
paper. Ref. 1. The assumptions made in that paper regarding flap geometry 
namely that the gap is sealed and that the flap leading edge coincides 
with the hinge line - are also assumed here. Summaries of the equations 
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and of the notation used by various authors will be found in textbooks, 
for example. Ref. 2. We prefer to use the notation developed by Wilts, 
Ref. 3, who clearly differentiates the physical origins of the terms. 

The airplane geometry and the associated coordinate systems are 
shown in Figure 1. The "reference axis" makes an angle A with the 
aircraft’s y a axis and lies at a distance * re £ behind the leading edge 
of the lifting surface. x re f m e& sure d perpendicular to the reference 
axis, x, z, 0, and 3 are structural degrees of freedom defined at a 
point on the reference axis. Aerodynamic forces and moments are applied 
to these same degrees of freedom. The flap hinge is parallel to the 
reference axis at a distance from the leading edge. The angular 
motion of the flap relative to the non-rotating part of the airfoil is 6. 
The positive direction of rotation is shown in Figure 1. The chord, c, 
and the spanwise dimension of the strip, Ar, are measured perpendicular 
and parallel to the reference axis, respectively. 

The equations for the aerodynamic pitching moment. Mg, lift, P^, 
and hinge moment, M^, all measured in the positive coordinate directions, 
are 
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The notation used in Equation 1 is defined in Table 1. The velocity 

V = V cos A (2) 

is the component of free-stream velocity perpendicular to the reference 
axis. The first term within brackets is the contribution of circulation 
and the remaining terms are respectively the non circulatory stiffness, 
damping and apparent mass terms. The coefficient matrices in the appa- 
rent damping and mass terms are 



The angle of attack due to motion at the rotation center, a » is 
the following function of the motions of the strip 

<W " 6 ♦ tan A 3 ♦ - i Cz ref - (C - x p - x ref ) 0 - cw 2 6) (5) 
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The angle of attack due to gust is 


a 

gust 


U(t-T) 

V 


( 6 ) 


where U(t-T) is the vertical gust velocity delayed by a time, x, which 
is specified independently for each strip. The coefficients w^ 
i= 1, 2, — , 14 are functions of. the hinge location, x^/c. They are 
closely related to Kussner's $ functions. Ref. 4. Formulas for the 
coefficients are listed in Table 2. The user will have the option in 
SADSAM to specify his own values of the coefficients. It is possible 
that some of the coefficients may be available from test data. For 
example, w^, w^, w^ and can be obtained from static tests in a wind 
tunnel. 

C(s) and ip(s), where s = cp/2V, are respectively the Theodorsen and 
Kiissner functions. They are approximated in SADSAM by the following bi- 
quadratic polynomial fractions 


(1 + d 1 s)(l + d 3 s) 
C(S) = (1 + d 2 s)(l “ d 4 s) 

1 + d^s 

= (1 + d 6 s)(l + d ? s) 


(7) 

( 8 ) 


The user has the option to select his own values of the coefficients 
dj to d^. The following default values have been shown to be good 
approximations for the two-dimensional incompressible case. Ref. 3. 


= 10.61 

d 3 = 1.774 

d 5 = 6.608 

dy = .834 

= 13.51 

d 4 = 2.745 

d 6 = 10.0 
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Chordwise Motions and Aerodynamic Forces 


Chordwise motions and aerodynamic forces are important for the aero- 
elastic analysis of fixed wing aircraft when low frequency stability and 
control characteristics are a consideration. They are important for 
rotary wing aircraft because the lower chordwise vibration modes occur 
in the same frequency range as the lower flapwise bending modes. This 
situation can also occur for some types of fixed wings. 

The coordinate geometry that will be used is shown in Figure 1. The 
x-direction is perpendicular to the reference axis. The following sketch 
shows the assumed flow condition in the x-z plane. 

z 


V = V - x 
f 

In the sketch, V is the component of the airplane's velocity 
relative to the free stream which is normal to the reference axis, and 
v is the induced velocity. For small angles, it may be assumed that the 
chordwise force on the strip is 

F x = D + L(J> (9) 

where D is the profile drag, 

D = ipV|cAr-C D , (10) 
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L is the lift. 


L = £pV*cAr*C L , 


(ID 


and <j> is the in-flow angle, 


(j) = 


z + v 


V f + a ind 


( 12 ) 


According to elementary momentum theory. Ref. 5, the induced angle of 
attack 

“ind " C L /lrA (13) 

where A is the aspect ratio of the surface. The induced angle of attack 
reduces the angle between the flow and the chord plane, so that, assuming a 
linear relationship between lift and angle of attack. 


C L = C Ha^ a ' a ind^ * 


(H) 


where C? is the lift curve slope for two-dimensional flow, 

X>(x 

of cu ^ between Equations 13 and 14 gives 


Elimination 


C L 


J £a 


ia 


1 + 


ia 

ttA 


(15) 


so that the lift curve slope including aspect ratio correction is 


J £a 


‘ia 


1 + 


tiA 


(16) 


Also, from Equation 14, 


, Sia\ _ C £a 
a ind "| » I nA 

Hot. 


( 17 ) 
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We are now prepared to write expressions for the forces in the z and 
x directions. Assuming <{) small and ignoring the small contribution of 
drag in the z direction 

F^ = L = ipV^cArC^ a *a (18) 

The forward velocity includes the dynamic perturbation x, i.e., 

V f = V - x (19) 

while the angle of attack, which includes both a steady part 0 q and a 
dynamic part, is 


“ ■ 0 o + C ^ a rot + ^“gust 


( 20 ) 


where a , a , Cfs) and i|/(s) are obtained from Equations 5, 6, 7 and 
rot gust T 

8, respectively. Note that the steady part of ^ is accounted for in 

the calculation of C„ and should not, therefore, be subtracted from 6 . 

ia U 

In the expressions for and the velocity V should be replaced 

by V^, but the effect of the substitution is second order small and will 
be neglected. If the lift is separated into a steady part and a 

i 

dynamic perturbation 6L, then Equations 18 to 20 give 


L 0 = iP V 2 cArC ta e 0 


( 21 ) 


and 


5F z = SL = ipV 2 cAr Cj, Jc (s ) c* rot + 'l'(s)a gust - 20 Q M + (nonlinear terms) 

( 22 ) 


The nonlinear terms will be neglected. The chordwise force given 
by Equation 9 can similarly be separated into a steady part and a dynamic 
perturbation. The dynamic perturbation is 
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6F =6 D + L n 6<j> + <ML 
x 0 Y T 0 


( 23 ) 


Lq and 6L are given by Equations 21 and 22. In addition. 


6D = - ±pj/ z cArC n *-#- 


a s to- 
^0 ttA 0 


t * 4r (C (»}«_* + 


Substitution into Equation 23 gives 


2C o 6n 
0 


6F x = ipV cAr — — (C(s)a rot ♦ *(s)a ) 


£ / < c *a e o )2 \i 

+ c 0 e n £ - 2 c n + r - — * 

Jla 0 V D ttA / y 


Equations 22 and 27 are used to generate perturbation aerodynamic forces 
for SADSAM. The perturbation pitching moment and hinge moment are obtained 
by multiplying 6F z by x re f x p and cw^ respectively, as in Equation 1. 

One further refinement is that the structural x and z axes may be 
rotated out of the aircraft's x a »y a plane through an angle y about the 
pitch axis. Thus, the forces applied to the structure are related to the 
forces computed in Equations 22 and 27 by 


cos y sin Y 1 (<$F 


■sin Y cos Y ( SF 
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while the velocities appearing in the equations are related to the 


structural velocities by 



— 




( 1 \ 

(M. 

cos 

Y 

-sin 

Y 

: 

Ui 

sin 

L. 

Y 

cos 

Y 

Uj 


( 29 ) 


The additional parameters required to describe the effects treated 
in this section are: 


0q, the steady angle of incidence of the lifting surface, relative 
to the free stream, which may be different for each strip, 

Cp, the drag coefficient, 

A, the aspect ratio, 

y, the rotation angle of the structural x and z axes, which may 
be different for each strip. 


Wing-Tail Interaction 


A lifting surface will induce an angle of attack on a lifting 
surface in its wake. This effect is important for the analysis of the 
response to gusts. The induced angle of attack can be expressed in the 
form 

where e is the induced flow angle (positive* downward) at the tail, and 
t 

C^(t-x) is the lift coefficient at some representative station on the 
wing, delayed by the time increment required for the wake to be convected 
from the wing to the tail. Hie dynamic part of will be included in 
the expression for angle of attack. Equation 5, for all affected tail strips 
but the steady part will Be ignored. (The user can include the steady part 
in 0 Q , if he wishes;.) ' 
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GUST ANALYSIS 


A gust analysis was made of a hypothetical small aircraft to demonstrate 
the gust analysis capability of SAD SAM. The definition of the aircraft con- 
figuration and the analyses were accomplished cooperatively with the con- 
tract monitor. 

Appendix A contains a definition of the aircraft configuration and a 
summart of some of the results of the analysis. 



I 

>-» 

t 



Table 1. Notation 


c chord 

s cp/2V, complex form of reduced frequency 

t time 


w 1 ,w 2> . 

x 


X 

p 

x ref 


♦ • a see Table 2 

14 \ 

chordwise coordinate, see Figure 1 for positive direction 
distance from leading edge to hinge line 
distance from leading edge to center of pressure 
distance from leading edge to reference axis 
vertical coordinate, see Figure 1 for positive direction 




C(s) 


D 

F 


x 

F 


z 


L 

L 

M 

M 


0 

6 

e 


U(t-T) 


V 

V 


aspect ratio 
drag coefficient 
lift coefficient 
lift curve slope 

lift curve slope for two-dimensional flow 
lift coefficient at some point on the wing 
Theodorsen^ function 
drag 

chordwise aerodynamic force 
vertical aerodynamic force 
lift 

steady part of lift 
hinge moment 
pitching moment 

vertical gust velocity at a particular strip 
free stream velocity 

component of free-stream velocity normal to the reference axis 
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Table 1. Notation (Continued) 


V, 



a 


a 


gust 

a ind 


a. . 

it 


a 


rot 


8 

Y 

6 

60 

<$F 

2 

6L 

6$ 


e 

e o 

V 

P 

T 

<J> 

*o 

^(s) 


angle of attack 

angle of attack due to gust 

induced angle of attack 

induced angle of attack at the tail due to lift on the wing 

angle of attack at the rotation center 

spanwise slope, see Figure 1 for positive direction 

rotation angle of the structural axes relative to the aerodynamic axes 

relative flap angle, see Figure 1 for positive direction 

dynamic part of drag 

dynamic part of F z 

dynamic part of lift 

dynamic part of <f> 

induced flow angle at the tail 

pitch angle, see Figure 1 for positive direction 
steady angle of incidence relative to the free stream 
induced velocity 
air density 

time delay due to convection 
inflow angle, see Equation 12 
steady part of 
Kiissner's function 


Ar spanwise width of strip 

A sweep-back angle 
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Table 2. W Functions 


Note: <p 


-£-) 


W 1 

= 

“ 

I 

IT 

0 + sin 0j 

W 2 

= 

1 

"4tt 

sin 0 (2 - cos 0) + 0 ( 1 - 2 cos 0) 

W 3 


1 

"5tt 

sin 0 (2 + cos 0) -0(1 +2 cos 0)J 

w 4 

- 

4 

IT 

sin 0(1+ cos 0)j 

W 5 

= 

ir 

0 + | (sin 0) (2 + cos 0)(1 - 2 cos 0)J 

W 6 

S 

l 

ff 

0-1 (pin 0) (2 - cos 0) (1 + 2 cos 0)1 

w 7 


4 

sin 0 (1 + cos 0) (0 - sin 0) 



TT 

2 ^ 

W 8 

= 

1 

”7 

IT 

|^sin 0 cos 0-0^ 

r l 

w 9 

= 

j, 

ir 

- sin 0 cos 0 + 0 

W 10 

- 

1 

77 

|^- 0 cos 0 + -j sin 0 (2 + cos 2 0)J 

W 11 

= 

1 

7¥ 

0’ cos 0 + i sin 0(1 + 2 cos 2 0)j 

W 12 

= 

1 

j^0 (A + cos 2 0) - ^(7 + 2 cos 2 0) sin 0 cos 0j 

*13 

s 

-L^f- 0 cos 0 + sin 0^0 - cos 0 sin 0j 

. W 14 

= 

_lF0 Z (I + cos 2 0) - ^0 sin 0 cos 0(7 + 2 cos 
4ff l _ -.1 


+ i sin 2 0 (4 + 5 cos 2 0) 


- 17 - 



REFERENCES 


1. Theodorsen, T. , "General Theory of . Aerodynamic Instability and the 
Mechanism of Flutter," NACA TR-496, 1934. 

2. Scanlan, R. H. , and Rosenbaum, R. , Introduction to the Study of 
Aircraft Vibration and Flutter, The MacMillan Co., N.Y., 1951. 

3. Wilts, C. H. , "Aerodynamic Forces in Analog Computation," California 
Institute of Technology, Computing Center, Technical Report No. 116, 
September, 1959. 

4. Kussner, H. G., "Zusammen fas sender Bericht iiber instationaren 
Auftrieb von Flugeln," Luftfahrt-Forsch. , 13, 412-424, 1936. 

5. Perkins, C. D., and Hage, R. E., Airplane Performance Stability 
and Control, John Wiley § Sons, Inc., N.Y., 1954. 


- 18 - 



APPENDIX A 

Aircraft Configuration Used in Sample Gust Analysis 


Wing : 



All Wing Mass CG 19" Aft of L.E. 


A-l 




Horizontal Tail: 





Stiffness : 

El = GJ = 3.07 x 10 7 LB. IN 2 

Mass : 


Sta 

In. 

Mass 

lb.sec 2 /in. 

Pitch 

Inertia 

lb.sec 2 /in 

0 

.0095 

.608 

33 

.0190 

1.216 

66 

.0190 

1.216 

99 

.0207 

.608 


Horizontal tail rotation about the 1/4 chord is restrained by a spring GO 
connected to the fuselage slope at Sta. 17, where: 

k = 1.18 x 10^ in. lb/rad 
A- 2 



1/2 Fuselage : 




2 

Mass(lb.sec / in) 

Stiffness 

? 

Sta. 


Empty 

Load 

El (lb. in ) 

Linear between Points 

- - 

-16 


.7617 

.0285 

4.62 x 10 8 

17 


.3135 

.1632 

4.62 x 10 8 
5 x 10 9 

47 

Wing 
& chord 

.3277 

.4741 

9 

5 x 10 

90 


.1736 

.3446 

5 x 10 9 q 
3.35 x 10 

130 


.0531 

0 

1.655 x 10 9 

175 


.0311 

0 

5.75 x 10 8 

220 

Horiz.Tail 
£ chord 

.0453 

0 

■ 

1.00 x 10 8 


SAS 

Horizontal Tail Rotation (6^) due to Pitch (6) 
« t (deg) = 

s 

S = 3 

{ t k e s 

T = i;T67s^iH;onT^T T (rad/rad;i ' 


0 = fuselage slope at Sta. 47 








Wing Rotation (^ v; ) due to Vertical Acceleration (z) 

<5 (deg) = k.. — \( z (in/sec 2 ) 

w z \s+0. l/\s+30 / , 

4 

5 

S = 3 

6 . 006464k.. s 3 

T = (3 . 33s+l )”(. 01 11 s+1) ( rad / in ) 

z - the vertical displacement of the fuselage at Sta. 47. 
The ’'basic” values of constants are: 


kg = .262 

k.. = 6.477 x 10 
z 


Stability Analysis Summary 


Case 

Run 

V 

SAS 


Stiffness 

Description 

knots 

k.. 

z 

k e 

factor 

5 





xl 

Normal Modes 

5 


100 

0 

0 

xl 

Flutter Analysis 

5 


200 

0 

0 

xl 

n ft 

5 


300 

0 

0 

xl 

n ft 

17 


170.7 

0 

0 

xl 

Stability of 

15 


170.7 

0 -4 

.262 

xl 

flexible airplane 

15 

2 

170.7 

6.477x10 l 

.262 

xl 

varying k" 

18 


170.7 

6.477x10 

.262 

xl 


18 

mm 

170.7 

.06477 

.262 

xl 


20 

If 

170.7 

6.477x10"!! 

.262 

xlOOO 

Stability of 

20 


170.7 

6.477x10 6 

.262 

xlOOO 

stiff airplane 

20 

n 

170.7 

.06477 

.262 

xlOOO 

varying k” 

16 


170,7 

-6. 477xl0 -4 

.262 

xl 

— 

Stability with k” 
sign inversion 


A- 4 












Frequency Response Gust Analysis 


The gust spectrum used is given by 

4>(a>) = (ft/sec) 2 /rad/ft 

4x10 +(,,j/U) 2 

The gust analysis was performed for an airplane velocity (U) of 
170.7 knots. The table of gust velocity (V ) versions of frequency (f) 

o 

used in the analysis was calculated from 

V (f) = ' - 9 - 6 - 5 - 1 - 1 — knots / &T 

g 8 . 4224x10" +f 2 

The printed or plotted output from SADSAM will bear dimensions of 

output units/ vhz 

For example accelerations, scaled in SADSAM to units of g, will be in 
units of 

g's/A\z . 


The mean square value, calculated by SADSAM for any output quantity 
4>, is given by 


f 



The mean square value of an acceleratiort output request scaled to 

2 

g's then has units of g . 



Run Description 

Case 21 Run 1 Flexible airplane, with SAS 

Case 21 Run 2 Flexible airplane, without SAS 

Case 21 Run 3 Airplane lOOOx stiffer, with SAS 

Case 21 Run 4 Airplane lOOOx stiffer, without SAS 

Output Summary 



Case 21 


Run No. 


Flex/Stiff 


SAS 


Mean Square Values 


2 

A22 z near airplane CG (g) 

A56 z at root of horiz. tail (g) 

2 

A2 z at wing tip (g) 

2 

A71 z at tip of horiz tail (g) 

•• 2 
A24 0 at root of wing (rad/in) 

A58 0 at root of horiz tail (rad/in) ‘ 

E4 , wing root bending moment (in lb) 

E28 horiz tail root bending moment 
(in lb) 2 

El 8 fuselage bending moment 
aft of wing (in lb) 2 

E20 fuselage bending moment 2 
forward of tail (in lb) 

2 

E14 wing root torsion (in lb) 

2 

E34 tail root torsion (in lb) 

2 

F7 gust load on wing cell (lb) 

2 

F66 gust load on tail cell (lb) 



Flex 


Yes 


. 3430E4 


.08063 

.05193 

.05691 

.7526 1 

.3535 

.3819 

1.3745 . 

.05152 

.05647 

1.2087 ! 

.3530 ; 

.3814 

.1136E-4| 

5.844E-6 

6.056 E-6 

1.986 E-41 

. 1314E-4 

6.059 E-4 

4.893 E8 

3.481 E8 

3.633 E8 

1,950 E6 

1.149 E6 

.9775 E6 

7.902 E6 

9.858 E6 

7.327 E6 

2.524 E6 

2 . 179 E6 

1.693 E6 

. 4692E6 

j . 3154E6 

. 3781E6 

. 2318E4 

| .1184E4 

778.3 

6.322 E4 

6.322 E4 

6.322 E4 

. 7351E4 

. 7351E4 

.7351E4 










Case 21 Run 1 


RMS Values 





A22 

.2751 

.2840 

.2279 

.2386 

A56 

.8963 

.8675 

.5946 

.6180 

A2 

1.1818 

1.1724 

.2270 

.2376 

A71 

1.0837 

1.0994 

.5941 

.6176 

, 

A24 

3.544 E-3 

3.370 E-3 

2.417 E-3 

2.461 E-3 

A58 

15.69 E-3 

14.09 E-3 

3.625 E-3 

2.462 E-3 

E4 

21431 

22120 

18657 

19060 

E28 

1602 

1396 

1072 

989 

E18 

3407 

2811 

3140 

2707 

E20 

1924 

1589 

1476 

1301 

E14 

641.9 

685.0 

561.6 

614.9 

E34 

58.57 

48.15 

34.41 

27.90 

F7 

251.4 

251.4 

251.4 

251.4 

F66 

85.74 

85.74 

85.74 

| 85.74 


RMS Ratio 


with S AS/with out SAS 

FLEX/ STIFF 

FLEX 

STIFF 

SAS 

No SAS 

.9687 

.9552 

1.207 

1 * 190 

1.0332 

.9621 

1.507 

1.404 

1.0080 

.9554 

5.206 

4.934 

.9857 

.9619 

1.824 

1.780 

1.0516 

.9821 

1.466 

1.369 

1.1136 

1.4724 

4.328 

5.723 

.9689 

.9789 

1.149 

1.161 

1.1476 

1.0839 

1.494 

1.412 

1.2120 

1.1600 

1.085 

1.038 

1.2108 

1.1345 

1.304 

1.221 

.9371 

.9133 

1.143 

1.114 

1.2164 

1.2333 

1.702 

1.726 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 











Transient Gust Analysis 


Airplane velocity ■ 170.7 knots 
Sharp Edge Gust 


50 ft/sec 


Case 25 - no structural damping 

Run 1 flexible airplane with SAS 
Run 2 flexible airplane without SAS 
Run 3 lOOOx stiffer airplane with SAS 
Run 4 lOOOx stiffer airplane without SAS 

Case 24 - structural damping added by means of SIMDAMP card in SADSAM, where 
g = .01 at 20 hz. Run schedule same as Case 23 above. 

Transient Output 

Same as frequency response solutions plus: 

V22 vertical velocity near airplane CG (in/sec) 

D39 pitch slope displacement of fuselage near the airplane CG (rad) 



General Notes : 


Consistent plotting scales between runs were accomplished for the frequency 
response solution by repeating the entire solution and specifying the plot scales 
chosen from a prior solution. This same task could be, but was not accomplished 
for the transient solution. 

Case 24 (transient solution) was performed with a small amount of structural 
damping. Case 23 was run with no damping and resulted in generally ’noisy 
solutions and in run 1, a high frequency instability. This high frequency 
solution content is probably the result of the numerical integration destabilizing 
some marginally stable high frequency root of the system. It is recommended 
that some structural damping be included in all transient solutions. Case 23 is 
included here to demonstrate that even the small damping included in Case 24 
(g s .01 at 20 cps) can have a profound effect. 

It should be noted that the gust loads on the horizontal tail are delayed 
.05 seconds after the gust loads on the wing. The effect of this delay can be 
seen in most of the plotted output and is most evident in the plots of F66, 
the gust load on a tail strip. In the frequency response solution the phase 
plot of F66 shows approximately a linearly increasing phase angle with fre- 
quency. The transient solution plot shows a .05 second delay before the onset 
of the load at F66. 

The transient solution converges to a constant angle of attack, equal and 
opposite to the gust angle (.1734 rad). 


Case 24 t = 1.0 seconds 


Run 

V22 
in/ sec 

D39 

rad 

V22 

knots 

n/v 

V22/170. 7 
rad 

a 

0 - n/v 

1 

228.92 

-.1070 

11.295 

.0662 

-.1732 

2 

122.87 

-.1380 

6.062 

.0355 

-.1735 

3 

231.80 

-.1062 

11.437 

.0670 

-.1732 

4 

128.78 

-.1363 

6.354 

.0372 

-.1735 
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